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ABSTRACT 


Investigation of urban effects on the power spectrum 
of turbulence is of interest in estimating the rate of 
digfuston of pollution inthe; city, and also in) determining 
the applicability of theories concerning atmospheric 
turbulence in a rough, non-homogeneous environment. In 
this study, turbulence observations from an urban site and 
a rural site at Edmonton, Alberta, were compared. Periods 
with light winds and stable lapse rates were chosen, because 
the rate of diffusion is more likely to be important under 
these conditions. The observations analyzed in this study 
were obtained on six days in September, 1970. 

The Pe reese renee were made using 
propellor-and-bivane anemometers (Vector-Vanes) mounted at a 
height of 2.7 m. Supporting lapse rate measurements were 
also obtained. The analogue signals from the Vector-Vanes 
were converted to digital form for analysis. Power spectra 
were computed by averaging short, modified periodograms. 
Turbulent intensities were computed and corrected for 
attenuation at short wavelengths. All turbulence statistics 
were computed in terms of polar coordinates. These statistics 
were used as approximations to their Cartesian analogues. 

It was found that the turbulent intensity was more 
than three times greater at the urban site than at the 


rural site. 
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The urban wind speed, elevation, and azimuth spectra 
agreed well with the Kolmogorov -5/3 law at wavelengths 
shorter than 8 to 50 times the distance to the ground. 

Rural speed spectra agreed with the -5/3 law at longer 
wavelengths than did rural elevation spectra. It was 
concluded that isotropy was not a necessary condition for 
a -5/3 law to apply to turbulence spectra. 

The position of the peak in the rural elevation 
angle spectrum agreed with similarity theory and recent 
observations. The position of the peak in the urban 
elevation spectrum was displaced toward longer wavelengths. 
Two alternative hypotheses were proposed: either that 
Similarity theory was not obeyed close to the ground ina 
rough urban environment, or that polar and Cartesian spectra 
were not similar in shape when the turbulent intensity was 


high. 
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CHAPTER it 
INTRODUCTION 


Meteorologists and geographers have been interested 
In che terrect Of a city on ine climate for some time. 
Studies have been carried out on the urban heat island, 
urban effects on precipitation, and urban effects on the 
mean wind field. One aspect of a city's influence on its 
climate that has only recently received attention is its 
effect on the magnitude and scales of air turbulence, 

A knowledge of the characteristics of turbulence 
in the caty 7s usercud, in two major. ways: first, in: 
estimating the effects of wind loading on tall buildings 
(Davenport, 1967); and second, in estimating the rate of 
diffusion of pollutants. This study is directed toward 
providing some data applicable to the second of the above 
problems. An experimental comparison of urban and rural 
turbulence will be made under conditions when the rate of 
diffusion of pollutants would be expected to be important 
to the inhabitants of a city. That is, when diffusion 
occurs slowly; in stable air with light winds. 

The relationship between the statistics of 
homogeneous, stationary turbulence and the diffusion of 
particulate matter in such a field was treated theoretically 


by Taylor ine P92). He welated the rate of ditfusion ‘to the 
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Z 
intensity of turbulence and the Lagrangian auto-correlation 
coefficient. Because the auto-covariance and the power 
spectrum are Fourier transforms of one another, it is 
sufficient to specify the Lagrangian power spectrum in order 
to estimate the rate of diffusion. In practice, the rate 
of diffusion has been measured and used to estimate the 
Lagrangian power spectrum since the latter is very difficult 
to obtain by direct means (Lumley and Panofsky, 1964). 

Various attempts have been made to relate the 
Lagrangian power spectrum to the Eulerian power spectrum, 
which is easier to measure. Among these studies is that of 
Hay and Pasgquill (1957), who suggested that Eulerian and 
Lagrangian spectra were Similar except that the time scales 
were related by a constant factorsectsingsthissa hypothesis, 
the Eulerian time spectra to be presented in this study may 
be, related to > the, rate) of diffusion of pollutants. 

A previous study comparing urban and rural 
turbulence is that of Bowne, Ball, and Anderson (1968), at 
Fort Wayne, Indiana. They found that turbulent intensities 
were greater in the city than in the country and that there 
were some systematic differences in the spectra. 
Consequently, there may be some significant differences in 
the diffusion of pollutants in rural and urban environments. 
Also, analysis of these differences may shed some light on 
the range of applicability of similarity theory in predicting 


turbulence statistics in varying environments. 
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The intensity and scales of turbulence may be 
affected by gross meteorological conditions in ways that 
are not completely understood. In order to isolate the 
effect of the urban environment, concurrent observations 
of the three components of the wind were obtained at an 
urban and a rural site. The observations were taken 
between 0500 and 0700 Mountain Standard Time, when 
inversion conditions were likely to prevail. Fast-response 
propellor-and-bivane instruments ("Vector-Vanes") were 
mounted on 2.7-m portable masts at both sites. Supporting 
estimates of lapse rates were made using a tethered balloon 
and data available from thermometers at various levels ona 
downtown office building. From these data, urban and rural 
turbulent intensities and power spectra were calculated and 
compared. The results to be presented were obtained on six 
days in September, 1970. 

ChapterodIvcontains’ a detailed description: of the 
collection and recording of the data. Mean wind speed and 
lapse rate tables which are used in later analyses are 
presented in chapter III. In chapter IV, the method 
used to calculate power spectra is described and the urban 
and rural spectra are compared. Chapter V is concerned 


with the calculation and comparison of turbulent intensities. 
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CHAPTER: Wee 
COLLECTION. OF THE DATA 
2.1 The Observing Sites 


Zain TheyRuralrSite 

The rural site was located on the University of 
Alberta farm at Ellerslie, about 5 km south of the city 
limits. / [ts position jisAindicated von yrigure = 1/ The 
instruments were situated near the eastern edge of a field 
of mown hay about 10 cm high. A band of trees bounded the 
field on the west. These trees were approximately 15 m 
high and 400 m from the observing site. A group of farm 
buildings was about 100 m to the north. Trees also occurred 
about 200 m to the) southeast; and to the northeast. .The 
ground slopes downward to the northeast at less than l 
percent. The only significant topographic feature nearby 
is the Whitemud Creek gully, which may be seen on the aerial 
photograph “(Figure 2): /It runs in a northeasterly direction 
and liesmat a distance of MWikm-to the west..of the site. {the 
Site is felt to be typical of the cultivated parkland of 


Central Alberta. 


A.062, The Urban Site 
In 1970 the population of the city of Edmonton was 


435,000; the built-up area was about 16 km square. The 
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Figure 1 Map of Edmonton and vicinity 
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Figure 2 Aerial Photograph of the Rural Site 
(Government of Alberta) 


etic pativrsedd + 


(sind orft. to dgsupotodd Istaeé §& siupit 
(stisdiA to ag lll 


North Saskatchewan River valley, which is 50 m deep, 

bisects the city from southwest to northeast. In view of 
the high artificial roughness and inhomogeneity of the city, 
its topography is not thought to be a significant influence 
on the low-level turbulence. 

The observing site was on a concrete plaza north of 
101A Avenue between 100 and 99 Streets. This location was 
near the geometric centre of the city and on the eastern 
fringe of the most densely built-up area. The plaza was 
30 m by 75 m with the long axis running east-west. The 
plaza was bounded on the north by the six-story Edmonton 
Public Library. The instruments were placed 25 m south of 
the library, approximately centered between two 9-m square 
one-story buildings in the southwest and southeast corners 
of the plaza. There was a 10-story building 75 m to the 
west of the instrument site. The buildings to the south and 
east were of 1 to 4 stories and their height dropped off 
gradually to the east. Figure 3 is a sketch map of the 


urban site. 


2.2 Turbulence Measurements 

Turbulence measurements were made using two 
three-dimensional anemometers (Vector-Vanes) manufactured by 
Meteorology Research, Inc. The anemometers consisted of a 
bivane and a propeller mounted on a shaft 51 cm long. The 
wind speed was measured by a light-beam chopper attached to 


the propeller; elevation and azimuth angles of the vane were 
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9 
measured by potentiometers. The output was in analog form. 
The pertinent manufacturer's specifications for these 
instruments were: 

Propeller response distance L = 2 - 3 feet (0.6-0.9 m) 
Vane delay distance.-D = 2 - 3 feet 

Vane damping coefficient o = 0.4 - 0.7 

The anemometers were mounted on portable masts so 


that the vane pivot was 2.7 m high. 


2.)3° “HSeamation. OL Lapse Rates 


2.3.1 The Rural Site 

Temperature measurements were taken at approximately 
30 m intervals from 1 to 150 m above the ground, using a 
thermistor bead suspended from a tethered balloon... The 
system consisted of a mercury battery and Wheatstone bridge 
to measure the resistance of the bead by the null method. 
The time constant of the bead proved to be too fast for a 
reliable null to be obtained. The system was therefore 
recalibrated using the imbalance current at four standard 
resistance settings to indicate bead temperature. Simple 
linear regression was used to fit the calibration data; the 
standard error was 0.1 C and the slope of the calibration 
curve was approximately 2.1 C mene A pilot-balloon 
theodolite was used to measure the angle of the balloon 
tether from the vertical. The height of each measurement 


was calculated from the length of the tether and this angle. 


No allowance was made for sag in the tether. This sag 
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(maximum difference between the chord and the arc) was 
estimated to be less than 10 percent of the arc length. A 
10 percent sag would introduce a systematic underestimation 
of the absolute value of the lapse rate of approximately 
2epercent. 

Two flights were made during most runs with measure- 
ments taken on both ascent and descent. Measurements at 
each level were approximately symmetrical in time about the 
mean time for the ascent. If the temperature at each level 
changed slowly during a run an unbiased estimate of the mean 
would be obtained by this technique. In view of the 
probable error in each measurement of approximately 0.2 C, no 
attempt was made to measure changes in lapse rate with 
time or height. All the data for each run were fitted by 
a straight line using simple linear regression. The standard 


error of the slope averaged about 10 percent. 


2.5.3 ‘The Uxban Site 

Geoscience Research Associates Ltd., under contract 
to the Department of Health, Government of Alberta, have 
taken vertical temperature difference measurements from the 
CN Tower building at 104 Avenue and 100A Street since 1967. 
The CN Tower is about 400 m north of the urban observing 
Site. Data from the CN Tower, along with data from other 
towers across Canada, are published regularly (in the 
Meteorological Tower Bulletin) by the Atmospheric Environment 


Service of Canada. The data used here are one-hour averages 
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¥d 
of the temperature difference between instruments at the 
17 m and 112 m levels. The thermometers are shielded, 
aspirated resistance bulbs mounted on 2.6-m booms. The 
accuracy of temperature differences between the sensors is 
estimated to be +0.1 F (Hage and Longley, 1970). The 
representativeness of the measurements is open to doubt in 
view of the proximity of the sensors to the building. The 
degree of horizontal inhomogeneity in the city is also 


unknown. 
2.4 Recording of the Data 


2.4.1 Analogue Recording 

The analogue signals from the Vector-Vanes representing 
wind speed, azimuth, and elevation were recorded by Precision 
Instrument.analogue F. M.. Tape recorders at. 1 7/8 i.p.s. 
(inches per second). At the beginning of each run the 
Vector-Vane output was disconnected from the tape recorder 
and zero Signal was recorded for about 30 sec. Zero 
output voltage corresponded to 0° azimuth, -60° elevation, 
and zero speed. The Vector-Vane was then reconnected and a 
"full scale" signal was recorded for another 30 sec. 
Full scale output voltage corresponded to 356° azimuth 
(because the Vector-Vane had a dead spot between 356° and 
360°), +60° elevation, and a speed of 20 miles per hour. 
The full scale voltage of the urban Vector-Vane was set at 
lea Volts and: that of “the rural Vector-Vane was set at 


4.0 volts. On playback, the gain of the tape recorder was 
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adjusted so that both the urban and the rural. signal had a 


BanGewors.s .24evolts. 

At the end of each two-hour run the procedure was 
repeated. These calibration marks went through the stages 
of filtering and analogue-to-digital conversion along with 
the data, and provided the relationship between the digital 
Signals and the original quantities speed, elevation, and 


azimuth. 


2.4.2 Analogue to Digital Conversion 

The data were converted from analogue to digital 
form by a Redcor’ Model 720 converter. The output-of the 
converter waS in increments of 5 millivolts. The sampling 
rate was 200 sec; four channels were multiplexed (one 
blank) so that the effective sampling rate for each channel 
was 50 sec and the Nyguist frequency was 25 Hz. 
Corresponding digital values for different channels were 
thus not simultaneous; a phase shift linear with frequency 
was introduced. The phase shift at the Nyquist frequency 
was 45° between adjacent channels and 90° between 
non-adjacent channels. 

The analogue signal was played back through low-pass 
filters to avoid aliasing. Adjustable active filters with 
an attenuation rate of 24 decibels per octave were used. 
The cutoff of the filters was set at 18 Hz, or about 0.7 


times the Nyquist frequency. 
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The analogue tapes were played back twice through 
the Gigitizer; once: ati1i)//38"1.<p.s., the same speed at 
which they were recorded, and once at 30 i.p.s., which 
compressed the time scale by a factor of 16. For this 
second digitizing run the effective sampling rate was 3.125 
sec; the Nyquist frequency was 1.56 Hz, and the effective 
PUVeer Seuvont ewas .lel25 hz. siThe) first«digatizing: run widl 
be referred to as the high-frequency run, and the second 
as the Low-frequency onun. 

Because of the limited frequency response of the 
Vector-Vanes at high wavelengths, combined with the low 
wind speeds encountered during the trials, almost all of the 
information available from the analogue turbulence signal 
was contained in the low-frequency digitized data. The 
high-frequency digitized data proved useful in demonstrating 


that aliasing had been avoided in the low-frequency data. 
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CHAPTER bid. 


MEAN WIND SPEEDS AND LAPSE RATES 


DURING THE TRIALS 


Lapse rates during the six trials, estimated by 
the methods described in Chapter II, are given in Table l. 
On September 4 the electrical connection to the tethered 
balloon at Ellerslie became intermittent due to the strain 
on the cable. Only two measurements, at 1m and at 50m, 
were made before the breakdown. These measurements indicated 
a decrease in temperature of 0.1 C from bottom to top of the 
layer. The uncertainty is very large, so this value is 
omitted from the table. No data were available from the 


CN Tower for September 13. 


TABLE 1 


Lapse Rates during the Trials 


Date Lapse Rate (C/100 m) 
Rural (Balloon) Urban (CN Tower) 
Sept. 4 ue Tlie o 
5 -2.1 -0.2 
9 -3.2 -0.8 
3 -6.7 a 
7] -5.5 -0.8 
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It may be seen from the table that an inversion was 
present during five out of five trials at Ellerslie, 
excepting September 4. Inversions were present downtown on 
three days. The lapse rate downtown on September 20 was 
Slightly less than the dry adiabatic, and on September 4 it 
was Slightly greater than the dry adiabatic. The lapse rate 
downtown was more unstable than that at Ellerslie on the 
four days when concurrent observations were available. 

Mean winds given in Table 2 were obtained from a 
10-m cup-and-vane anemometer at Ellerslie and from the 
Atmospheric Environment Service observations at Edmonton 
Industrial Airport (XD), which is located in the northwest 
quadrant of the city. It may be seen that the more unstable 
urban lapse rates of SepLGHDES 4 and September 20 were 
associated with higher wind speeds than occurred during the 


other four’ trials. 


TABLE 2 


Mean Wind Speed at 10 m 


Date Rural (Ellerslie) Urban. (XD) 
Direction Speed Direction Speed 
(2 SEYrus.) (m ea) (° true) (m eaor) 
Sept. 4 S15 BA 300 Led 
5 150 CRLT 180 4.0 
9 200 Sil 190 3.6 
is 180 22 190 Zee 
aly 210 ea 210 Bil 
20 279.0 ORES | 280 4.9 
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On September 4 the wind speed at Ellerslie was 5.4 
m sec +, while on the other five days rural wind speeds 
were between 2.2 and 3.1m gon ®. The lapse rate observed 
downtown on September 4 was super-adiabatic. Since strong 
winds are generally accompanied by breakdown of the surface 
inversion layer, it is very likely that the rural lapse 
rate was much greater on September 4 than during the other 
trials. THis: conelusien: is supported by the limited data 
available from the tethered balloon on that date. The 
conclusion that no inversion was present on September 4 at 
Ellerslie will be used in chapters IV and V. 

The mean wind speed at 2.7 m was calculated from 
the low-frequency digitized Vector-Vane data. The mean 


longitudinal component u is related to the mean speed s by 


lag Seo St sae es eae (2015) 


Zar 2 
to second order, where on is the azimuth variance and ots 


is the elevation variance, both in radians. When the 
turbulence is intense (o,* and a.¢ are large), the mean 


wind speed is an overestimate of the mean longitudinal 


this relationship may be derived from the velocity 
transformation equations for polar to Cartesian coordinates 
(with the x-axis parallel to the mean wind): 
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Component. ‘Table S»contains values of s and u- calculated 
Grom (3.1) and the percentage by which Ss exceeds W. A 


significant difference between s and u occurred in the city. 


u = s cos e€ cos a 
V»=. S$ cos e sina 
Wi} = S stn e 
In the first-equation, let s,= s + s'.and expand the 


cosines in Taylor series, retaining only second order terms: 


u= (3 +s') (1 - e*) (1 - a’) 
eA Zz 


and average both sides of the equation, again discarding 
terms of higher than second order, which yields 


u=s(l - 


7 : 
‘es 


a 


ne 7 woo e 
co pa 


esi ot ase na oi omens | 


rie 


Wed 


. ages fa, <a? . 
a as 
ea: arene baie taal 
: wu -—-* oT ry Mn Pa oe 7 ye | 


i, Voge 8 apap pale, 
- 2 Ate 7 bat 
i ; a pu ears wey Ue 
wll 
s "v 
a ol 


& Son & 209 2 _ o. 
5 aie S09 B= we an 
; i. 


18 


8T 90°? Oo ay 0 LoL EG amt 07 
ee a as 0 Sigeac BC ac ET 
TC 6L°0 C0 “et 0 GO eC 90 ve oat 
on L620 1.0 “SE ae ie a 6 
Sib Le al LO ah 0 OT sc Or eee S 
97 eae 90° | Ce" OS°? 09°? yp °3des 
so a a 2 ae ee ee ee ee eee eee eee eS 
(jueo7Ed) (2 des uw) (6 00S Ww) (UDd070ed) (_ oes w) (@ Das 1) 
s/(n-s) S 3 s/(n-s) So & 8 
ueqin Teazny 37eqd 


i few te W pue.s 


€ ATA. 


Yvan 


CHAPTER IV 


POWER SPECTRUM ANALYSIS 


A le Methods and Limitations ‘of Analysis 


4.1.1 Averaged Modified Periodograms 

Power spectra of wind speed, vertical vane angle, 
and horizontal vane angle were calculated from the low- 
frequency digitized data using a method described by Welch 
(1967). The method involves dividing the record into 
sections and calculating a modified periodogram from each 
section. The power spectrum is obtained by averaging the 
periodograms. A short summary of the method is presented 
below, indicating the choices of such things as data windows, 
and length and number of sections, which were made in this 
application vel THe *notation is essentially that used by 
Welch. 

The record was divided into K sections, X,, k = 
seny Ke Bach, section-consisted of L consecutive points. 
In this study K was either 38 or 42, L was 1024, and the 
sections overlapped by one-half their length. 

Each section was multiplied by the Hanning data 
window: 
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and the discrete Fourier transform was taken: 


ee 
Gia Beek eR A) EWG)eexp (-3. 2njn/L) 
g=0 
oe oe ee a = 0 bow a2 (4,2) 


The frequency n is expressed in cycles per unit time, 
where the unit of time is the length of one section of 
thes record, pore 328 sect. 
The Fourier transform was performed using a Fast 
Fourier Transform subroutine available in the University 
of Alberta program library. The algorithm, which was 
written by G. Sande, is described in detail by Cochran et 
alter (L967) the length of the series, L, must be an integral 
power of 2. 
Modified periodograms were calculated from each 


Fourier transform: 


L 
where Vie ys el y Wiretass (4.4) 
5 


Division by U is necessary to correct the variance for the 
effect of pre-multiplication by the data window. The value 
of U for the Hanning data window is (approximately) 3/8. 


The factor of 2 arises from the definition of spectral 
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density in terms of positive frequencies only. The power 


spectrum was obtained by averaging the periodograms: 
K 
Sia) (S=n)k ET, (n) (4.5) 


The spectral window corresponding to the Hanning 
data window is given by 


_, Ll 2 
Q(n) = U L © Wj) exp(-i 2ajn/b (4.6) 
dl E 


When the area under the curve is normalized to unity, the 
half-width of this spectral window is 0.8 and the height 


Or the first. sidelobe ius SL0r times that of the main 


lobe. 
A measure of the stability of the estimates is 
Var [S(n)] oaks | K-1 
——_—————- = kK lat 2nd (K-53) 9 (3) (4.7) 
BE [s(n)] JEU 
where 
moo | arene ne ae 
o(j) = om W(k) W{k+ 3 (L-M) } vi W (k) (4.8) 
k=0 k=0 


and M is the amount of overlap (in this case M = L/2). 
For n = 0 and n=L/2 the variance is greater by a factor of 
2. The quantity o(j) expresses the degree of dependence 


between the sections of the record when the spectrum is flat. 


(2,8) 
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For the same length of record, K is almost twice 
as large when the sections overlap by one-half their length 
as when they do not overlap. The resolution is unaffected 
by overlapping, so this technique uses the available data 
more efficiently. Greater stability is achieved without 
loss of resolution. 

Because for the low-frequency runs K was 38 or 
greater, the standard deviation of each estimate is about 
14 percent of its expected value. 

Power spectra were also calculated from the data 
digitized ,ae?50 se¢ i? The method was essentially similar 
to that given above, except that the sections were 2048 
poeintse (4irsécyei6éng*anhdethey"did notfoverlapl=+A 
noise peak was found at 10 Hz, which may have been due to 
mechanical noise in the tape recorder. Because the noise 
peak was almost a decade above the filter frequency for the 
low frequency runs, it was attenuated by over 70 decibels. 
This amount of attenuation removes any possibility of 
aliasing on the low frequency runs. Aside from this 
demonstration the high-frequency spectra presented little 


extra information. The response of the Vector-Vanes did 
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not extend to frequencies much above 1 Hz, and the 
resolution was insufficient to indicate the shape of the 
spectra adequately at low frequencies. 

4.1.2 Response of the Vector-Vanes at Short 

Wavelengths 

At short wavelengths the power spectrum is 
eeeenliatedeby sthetresponsejoiwthe Vector-Vanes. Under 
favourable conditions the observed spectrum may be corrected 
for this attenuation, extending the useful frequency range 
of the analysis. 

The bivane which was used to measure the azimuth 
and elevation angles of the wind behaves as a second order 
system and the propeller behaves as a first order system 
(MacCready and Jex, 1964). The relevant energy transfer 


functions are 


—1 
T, (n) = E ston (2 1h 1, Ae | (4.10) 
for the first order system, where L is the wavelength 
corresponding to n and lL, is the response distance of the 


propeller, and 


for the second order system, where L, is the natural 
wavelength of the vane and op is the damping coefficient. The 


natural wavelength is related to the delay distance by 
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(MacCready and Jex, 1964) 
L = (620 = 2.469) D% (4.12) 


The transfer functions are plotted in Figure 4. 
The observed power spectrum S(n) is related to the 


"true" power spectrum S,(n) by 
Says =) ne Sent) (4.13) 


At a wavelength of 4m, T, (nm) varies from 0.50 for 
Drpete Oi eOsil 2G USO ae Olas lig = (0.8 mM, T. (n) varies from 1.56 
fovershoot) foriDr="0.8 mand o%= 074 te 0.52 for D =.0.8 m 
anc 6 = 0.7. The functions then grow rapidly smaller. At 
a wavelength of 2 m both transfer functions are of the 
order of ome. at wavelengths equal to the delay or response 


distance, T, (n) US ;abpoul. 2 ox os and Te (n) is an order 


2 
of magnitude smaller. 

The uncertainty in the transfer functions resulting 
from an inexact knowledge of the response and delay distances 
and the damping coefficient is seen to be large. For this 
reason, the spectrum at wavelengths shorter than 4 m was 
discarded. Wind tunnel tests to determine the characteristics 


of the particular vanes in use would be necessary before 


correction could be expected to be reliable. 
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apt. ONAN Cr zat LON Gr rTOor 

When an analogue signal is converted to digital 
form, noise” rs“ introducedk “This notse ts. a result*of the 
rounding off of the analogue signal to the nearest digital 
value. When the total variance is expressed in the units 
Or tne drgical scale, the CONtrloutLLon oL this” quantization 
moLse=is iy t2 Unrt™ (pendac and-=Pilersott, 906, D."281)%- The 
smallest variances measured (for horizontal vane angles 
at Ellerslie) were of the order of 10 digitization units; 
most variances were at least one order of magnitude larger 
than this. Under all circumstances the guantization error 
added less than 1 percent to the variance. It is therefore 
neglected 

The "eLrect’ or quantization Grror on measured power 
spectra is also of interest. If the quantization errors 
for successive points are independent, the spectrum of the 
quantization error will be flat. For the low-frequency runs 
the Nyquist frequency is 512, considering the time unit to 
be*the™=rengeh or” one 1recora;, o1r* 526 sec. “"Iheretrore the 
power spectral density, in digitization units, of the 
GUaNCLAaL ron, SuLOrr lS t/t zee ly Oo OLD Leon x tos 

The logarithmic spectra contained in the Appendix 
are expressed in these arbitrary units for comparison. It 
may be seen that at the highest non-dimensional frequency 
plotted (about 0.65, or a wavelength of 4 meters), S(n) was 
of order pga4 for most spectra; on September 13 and 17 at 


Ellerslie S(n) for horizontal vane angle was of order 100. 
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Théreforenquantizatronocerror)icontributedyinsignificantly..to 
S(n) except for the above mentioned horizontal spectra, 
which contained about 10 percent noise at high frequencies. 

At non-dimensional frequencies higher than 0.65, 
the measured power spectrum decreased rapidly with frequency 
because of attenuation by the Vector-Vanes. At f = 1, 
most spectral densities were of the order of 5 x ‘Oe ope 
less. About 30 percent of this spectral density may be 
attributede toe quantizationa errorzo) Correction) of -spectra-for 
the Vector-Vane attenuation, as discussed in the last 
section, would result in amplification of the quantization 
error as well and the power spectrum would be seriously 
overestimated. Reliable correction would require 
amplification of the signal before digitization to reduce 
the relative contribution of the quantization error. 

Other sources of electrical and mechanical noise 
were evident in the analogue signal. The electrial noise 
was estimated at 50 millivolts before low-pass filtering, 
whieh corresponded to: LOGdigi tization units... li sthe 
spectrum of this noise was flat up to the high-frequency 
limit of the tape recorder, which was about 300 Hz, the 
contribution of this noise to the total variance would. have 
been 0.5, in digitization units. The spectral density would 
have been about TOKE es This noise does not contribute 
Sigurt icantly, to the Localvwvariance and ar pects tthe power 
Spectrum only for rural azimuth spectra’ on september: 13 and 


17 at high frequencies. Correction of the spectra for 
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Vector-Vane response would, however, require the reduction 
of this noise because it would be amplified by a large 
factor at high frequencies. 

Consideration gofathe quantization error and other 
sources of noise, in addition to the effect of uncertainty 
in the- transfer coefficients, indicated that the corrections 
which were made to the spectra were unreliable. Consequently, 
the corrections were not applied in the final analysis. The 
rapid increase of the signal-to-noise ratio with frequency 
also indicated that success would be unlikely in correcting 
Vector-Vane spectra at wavelengths shorter than 2 m, even 


if the transfer functions were known precisely. 


ae2) Discussion voL Results 

The power spectrum estimates are presented as 
functions of a non-dimensional frequency defined by 
Gi= 072/U,) wheres! 1s frequency, gsic the height of the 
observations, and u is the mean longitudinal component of 
the wind. This definition was chosen because it represents 
the ratio of height to wavelength at all levels of turbulent 
intensity to which Taylor's hypothesis applies. Numerous 
authors (e.g., Panofsky and McCormick, 1954) defined a 
non-dimensional frequency by f = nz/V where V is the mean 
horizontal wind speed. These early studies were concerned 
mostly with the vertical component of turbulence, and 
estimates of u were probably not available. The definition 
of non-dimensional frequency used in this study has been used 


DY Other authors, (6.0.5 Munn, 1966, :p. -70).. 


| ia’ ie 
t ae an 
~ | fate pee 


testo bus xow8e | 
yinisdxs0nm Fo bland las 
anotivextoo st tede inasialiis 
.vidensopeenod olde bieway ers oon anon hi ‘ 
oft vJeteyiens fem? oft ot ie! ot 


vorteupstt dsiw oftax sa ton-o 


pnitosames ai yloxitau sd biuow ees 
neve \m & osdt todrore ncenapehiiid suitoege « 
. ,.-ybezioesq awornd ai iiasiens necsbiith eenliailiat oid 


e6 betneestq 216 iedenvamieces mesa os ad? 
oes 


, 
ant ~~ 


ia 7 

éit to siipred edt ef &s , Yoreupott = a peer hy Ne = 
ERA: 

to Jnenoqmoo Ismibutipmol dsem sis et eB é brs —_— - 


etrereirges ti sauBasd noaoits enw stoke tek exer poe 
jasiudwwt to efevel {fs 46 alate ot thoied “te: wot 
avoitsmuy -eoctegs aaa. 2 ‘s08 galt toss oF 


5 bentiteb (neer = L eo 1s wh oie Ronee A 


yd heniteS yonsrpsx2 Lenotememtb=aon 8 20 8 


reom ott oY sna Vga = at no 
7 — - y 7 : me rf . - rir = 

benteono> staw 2s yluse sesdT 5 | 
t een oR oi - ee yt 


bos ,eonsludies te Fhenoqnion Isoisn: we 
ina 7 "ach ve a aves eS. peeves ake . 


see _ 


Pas) 


Theaquanwities Vand u-are related by 


Ue = Vi Gletaacede) (414) 


to second, order. s:Consequently when ‘the turbulent intensity 
is low, the definitions coincide. When the turbulent 
intensity is high, nz/V is smaller than the ratio of height 
to wavelength. For example, when the lateral turbulent 
intensity is, 0.de) nz Vou sts per cen trsmalder! than* nz /u. 
Because atmospheric turbulence spectra extend over many 
decades the difference between the definitions is usually 
negligible. It should, however, be kept in mind when making 
comparisons with the results of other studies. 

Graphs of the power spectra are presented in the 
usual forms, with the logarithm of non-dimensional frequency 
as the abscissa. The ordinates are the logarithm of power 
spectral density S(n), and the product of frequency and 
spectral density n° S(n). The first ordinate is expressed 
in arbitrary units. The second ordinate is normalized by 

2 


u“*so that the area: under the curve.is equal to the turbulent 


intensity. 


4.2.1 Vertical Vane Angle Spectra 
The power spectrum of vertical vane angle (elevation) 
may be identified with the power spectrum of vertical 


velocity through the approximation 
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W = ue! (415) 


where w' is a fluctuation in vertical velocity and e' is a 
fluctuation in elevation. Thus the shape of the spectrum 
of vertical vane angle is approximately the same as that of 
the vertical velocity. The higher order terms which have 
been neglected in (4.15) are very complex, and it is not 
clear under what circumstances the spectra of vertical vane 
angle and vertical velocity may be of significantly different 
shape. 

The! graphs) of “no-iSdn)eversus log £,.(Figures.A-1l to 
A-6) were smoothed by eye and estimates were made of Ea! the 
non-dimensional frequency for which n + S(n) is a maximum. 


The values of i which were obtained are given in Table 4. 


TABLE 4 


fn for Vertical Vane Angle Spectra 


Date Rural Urban 

Sept. 4 Ois-2 Q...16 
5) 0.4 0.08 
9 a 0.06 
Ms 0.33 Oie1.0 
W7 0.4 = 
20 Oes Danl2 


average 2.6 O10 
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It was originally proposed by Gurvich (see, e.g., 
Lumley and Panofsky, 1964) that according to similarity 


theory, 


Neos oSs (mn) + “gab (nZAvy Rs) (4.16) 


where Ss. (n) is the power spectral density of vertical 
velocity and F is a universal function of non-dimensional 
frequency and Richardson number. Therefore fn should be a 
universal function of Ri only. Busch and Panofsky (1968) 
found that fn was 0.32 under unstable and near-neutral 
conditions and increased rapidly with increasing stability. 
Their results were based on many power spectra from different 
locations in the past decade. Stability was measured as 
the ratio z/L, where L is the Monin-Obukhov stability 
length (Monin and Obukhov, 1954). This ratio approaches 
zero near the ground, and is zero at.all heights in neutral 
air. Consequently fon should vary little with lapse rate 
close enough to the ground, and the value of fn under these 
conditions should be the value for near-neutral air. 
Smoothed and averaged spectra for Hanford, Washington, at 3 
and 6.1. mishow.dittile) variation: with, stability (Busch and 
Panofsky, 1968). The average value of fn at Ellerslie (0.3) 
agrees well with the results of Busch and Panofsky for 
near-neutral air. The Hanford results indicate that fn at 
2.7 m should have the near-neutral value. Consequently, it 
is concluded that the value of fn obtained at Ellerslie is 


in agreement with the similarity theory. 
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The average value of =o at the urban site was 0.1, 
under conditions when the lapse rate between 50 and 150 m 
ranged from unstable to stable. This value is not consistent 
with the universal function proposed by Busch and Panofsky 
(1908) ,“byy which fa would never be less than 0.32 under any 
Stability. It is therefore suggested that, in the rough 
urban environment, vertical velocity spectra may not obey 
Similarity theory. This conclusion is not surprising since 
two assumptions made in the similarity theory are that the 
height is large in comparison to the roughness height, and 
that the terrain is homogeneous (Monin and Obukhov, 1954). 

If the earlier definition of f = nz/V were used, 
values of fn obtained at Ellerslie would be unchanged, 

Since the turbulent intensity was always low. Values of 

fn obtained at the urban site would be lower, since u was 
estimated to be an average of 14 percent less than V 
downtown. The average fn found using this definition was 
0.086. The conclusion that the observed values of fn 
downtown are inconsistent with similarity theory is therefore 
strengthened rather than weakened by the choice of the 
definiti on™ ofs & 

The possibility also exists that the relationship 
between the spectra of vertical velocity and vertical vane 
angle broke down. The relationship appeared valid for the 
Ellerslie spectra, but turbulent intensities downtown were 
much higher and the terrain was inhomogeneous. Regretfully 


tae question must be left open for the present. 
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Bowne, Ball, and Anderson (1968; also Bowne and Ball, 
1970), in their comparison of rural and urban turbulence 
at Fort Wayne, Indiana, concluded that vertical velocity 
spectra showed some shift of energy to shorter wavelengths 
in the city. Their measurements were taken at greater 
heights (about 15 and 60 m) when the rural stability was 
generally near-neutral. They ascribe the shift toward 
shorter wavelengths to an increase in mechanical turbulence. 

The results of the present study indicate a shift of 
energy toward longer wavelengths in the city, at a height 
of 2.7 m. Since the shift occurred on days when inversions 
were present at both locations, and since the observations 
were made close to the ground, stability would seem to be a 
minorafactors | tifsthisrshiftrmayiakisovbewascribed to 
mechanical turbulence, the scale of this turbulence is of 
interest. The horizontal scale of the mechanical turbulence 
is about 10z in Edmonton at 2.7 m: A typical spectrum for 
Fort Wayne at 53 m (Bowne, Ball and Anderson, 1968, p. 65) 
indicates maximum energy at 1.25z. A consequence of 
Similarity theory, which is obeyed reasonably well in the 
country, is that the scale of turbulence is proportional to 
height. Combination of the results from this study with 
those of Bowne, Ball, and Anderson indicates that this is 
not the case in the city. Rather, it is indicated that the 
scale is larger near the ground and smaller aloft than in 
the country; and therefore, increases more slowly with 


height. 
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4.2.2 Horizontal Vane Angle Spectra 
To gairsk- ondergtad latenaltvelocityef luctuationnvy« 


and an azimuth fluctuation a' are related by: 


Vinsetru av (4x7 ) 


Therefore, the spectrum of lateral velocity and the 
azimuth spectrum are of approximately the same shape. The 
approximation may be expected to be better than the vertical 
velocity - elevation approximation because wind speed and 
lateral velocity are not correlated in the same manner as 
wind speed and vertical velocity (Lumley and Panofsky, 1964, 
p. 175). Azimuth spectra will therefore be interpreted as 
lateral velocity spectra (normalized by ait 

The time series of horizontal vane angle contained 
discontinuities where the angle jumped from 0° to 360°. 
Sections of the records containing these discontinuities 
were discarded for analysis. During all trials except 
September 20, the urban horizontal vane angle record 
contained so many discontinuities that only one or two 
sections, about 400 sec long, were available for analysis. 
Consequently, only the low-resolution analysis made from the 
high-frequency run had adequate stability. This analysis 
did not extend to low enough frequencies to determine fn 
As a result only one value of fn is available for the urban 


site and no conclusive comparison of tn can be made. 
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On September 20, tn was 0.025 downtown and 0.063 at 
Ellerslie, which indicated a shift of energy toward longer 
wavelengths in the city (Figure A-13). 

Bowne, Ball, and Anderson (1968) found that lateral 
velocity spectra were roughly similar in shape in urban and 
rural environments at Fort Wayne. They suggest that there 
was some tendency toward more energy at shorter wavelengths 
iimeneacicy . 

The shifts observed by Bowne, Ball, and Anderson and 
during the present study are opposite in sense, as was the 
case for vertical velocity spectra. There is insufficient 
evidence to conclude that the lateral scale of turbulence 
behaves similarly to the vertical scale in the city (that 
is, increasing more slowly with height). The evidence 
available does raise this possibility, however, and 
further investigation would be of interest. 

An interesting phenomenon took place at Ellerslie 
during some observations. On September 5, 13, and 17 the 
azimuth power spectrum tended to increase towards lower 
frequencies after passing through a minimum in the 
neighbourhood of f = 0.01 (Figures A-11 and A-12). 
Examination of variances calculated from means over different 
periods supported this tendency. Table 5 shows the 
percentage increase in variance that resulted when the 
means were calculated over 656 sec rather than 328 


sec. 
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TABLE 5 


Percentage of Rural Azimuth Variance between 
OR OGLS32H2Z sand 0.00506. Hz 


Date 


Percentage 
Sept. 4 S 
S 32 
sige: 10 
i, 42 
20 26 


These values are roughly the amount of variance 


3 


attributable to frequencies between 1/656 Hz (1.53 x 10 ~ Hz) 


& 


and 1/328 Hz (3.06 x 10 ~ Hz). These frequencies correspond 


to values of f£ of 2x zone and 4 x loge and wavelengths of 
1080 and 540 m respectively, when the wind speed is 2 m sec >. 
These estimates were calculated from about 10 sections 656 
sec long. The variance between these frequencies may be 
interpreted as a single spectrum estimate S(n) at some 
intermediate frequency. The stability of the estimate may 


be obtained from equation (4.7). Because the sections 


did not overlap, 


Var [S(n)] 
> eae (4.18) 
me = bS-(1): | 10 
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Thus the standard deviation of each estimate is about 30 
percent of its expected value. In this light, there appears 
to be a significant difference between the low-frequency 
spectrum on September 4 and that during the other four trials. 
An inversion was present during the four trials when the 
phenomenon was observed, and the wind speed for the four 
trials was about 2m sect. It was shown earlier that the 
wind speed on September 4 was higher than that during the 
other trials, and that the inversion present on the other 
days was very likely absent. 

It may be concluded that significant power was 
contained in fluctuations with wavelengths between 540 and 
1080 m under stable, light wind conditions. 

This phenomenon is most, likely the same one that is 
mentioned by Lumley and Panofsky (1964, p. 178). They 
speak of "slow-moving features that may have horizontal 
dimensions of several 100 (sic) meters to a few kilometers", 
of unknown cause, which have been observed in stable air 


with light winds. 


4.2.3 Wind Speed Spectra 

Because fluctuations in the wind speed and in the 
longitudinal component are equal to first order, the shapes 
of their spectra should be the same. As in the case of 
vertical spectra, the conditions under which the first-order 
assumptron: breaks down areynot clear. » Consequently the 


approximation should be made with caution. 
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Comparison of urban and rural speed spectra, which 
are. given,in the Appendix, (Figures A-5 to A-10), does not 
reveal any obvious,difference: suchas; that:.foundfor 
vertical spectra. ,»~Valvueseot fn were obtained from the speed 


spectra by visual inspection, and are given in Table 6. 


TABLE 6 


fn for Speed Spectra 


Date Rural Urban 
Sept. 4 Ore Onl: 05.03 
5 0.06 0-025 
9 : 0.04 
NS double peak 0.06 
Le no well-defined peak 
20 Or 0cbG On 025 


On September 4 and September 20 there was a small 
shift of energy toward higher frequencies downtown, but the 
level. of Signiticance: was, ow. On September 5 the shitt was 
in the opposite direction. On September 13 there appeared 
to be two peaks in the spectrum at Ellerslie. Consequently 
no systematic. differences between rural and urban tn values 


are evident. 
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In the same manner as in the analysis of the 
horizontal spectra, the difference between variances 
calculated from 328-sec means and 656-sec means was 
examined. Again, the difference in variance may be 
attributed to wavelengths between 540 and 1080 m when u is 


Tn estes These differences are shown in Table 7. 


TABLE 7 


Percentage of Speed Variance between 
On00RS3FHzhand{ Ot00B06eHz 


Date Rural Urban 

Sept. 4 6 0 
is) One ibe) 
9 oF Se, 
13 4k, 3 sal 
Li Z> ais 
20 25 O97, 


Examination of the table shows that on September 17 
and September 20, the behavior of the rural wind speed 
spectrum was Similar to that of the lateral velocity 
Spectrum. That is, a large proportion of the variance was 
contained in eddies of wavelength between 540 and 1080 


vine A Similar phenomenon has been observed by other 
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investigators (Lumley and Panofsky, 1964, p. 184). 

The energy between 0.00153 Hz and 0.00306 Hz was 
low for all urban trials. The absence of the "meandering" 
observed at the rural site may have two causes. The first 
is the decrease in stability which was observed in passing 
from the rural to the urban environment. The second is 
concerned with the location of the urban site, which was 
actually between the major roughness elements (the surrounding 
buildings). The spacing of the buildings may place an upper 
limit on the wavelength of the disturbances to the flow 
between them. Further observations using refined techniques 
for lapse rate measurement would, hopefully, distinguish 
between these alternatives. 
4.2.4 Behaviour of the Spectra at High 

Frequencies 

It was first proposed by Kolmogorov (1941) that, 
when the Reynolds number is large, there should exist a 
locally isotropic inertial subrange of the spectrum in which 
no production or dissipation of energy takes place. The 
three-dimensional spectrum E(k) is then determined by the 
wave number k and the total dissipation e€ only. By dimensional 


reasoning it follows that 


ree ae eerev aS: (4.19) 


where A is a universal constant. Because of the requirement 


of local isotropy, all three one-dimensional spectra must 
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also be proportional to goes Also, because k = fz, S(f) 


7s The , lowerglimattofinthesinertial 


iswonopembionabrté, fi 
Ssubrange is expected to be at wavelengths short in 
comparison to the distance to the ground (Lumley and 
Panofsky, 1964) and short in comparison to the Monin-Obukhov 
scale:length.. Since the distance: from:’the ground: in’ the 
present study is expected to be small compared to the scale 
length, the first criterion will apply. 

Observations have indicated that one-dimensional 
spectra obey a -5/3 power law at wavelengths longer than 
predicted by the theory (Lumley and Panofsky, 1964). 
Evidence from the present study supports this conclusion. 

In fact, the short-wavelength limit of the power spectra 

in the Appendix is about twice the distance to the ground. 
Consequently any agreement with the -5/3 law that is present 
is at longer wavelengths than is explained by the theory 

of the inertial subrange. 

Logarithmic plots of spectral density against 
non-dimensional frequency are shown in the Appendix. The 
theoretical -5/3 law is indicated on each graph by a solid 
line. It may be seen that the logarithmic spectra generally 
agreed with a -5/3 power law at high frequencies. The fit 
was best for speed spectra, which exhibited small scatter 
over aidarge ranges of frequency. The rural-speed spectrum 
on September 13, however, showed an unexplained increase in 
slope above f = 0.1 (Figure A-21). The elevation spectra 


showed more scatter than the speed spectra (Figures A-14 to 
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A-23). The rural azimuth spectra agreed with a -5/3 


law on September 5, 17, and 20, but did not on September 4 
andst3 @thrguines A-249to A-26))); 

The urban azimuth spectrum on September 20 (Figure 
A-26) also fitted a -5/3 law well. The high-frequency, 
low-resolution urban azimuth spectra available for other 
trials fitted the -5/3 law well, but the resolution was 
insufficient to determine the wavelength at which the law 
no longer held. These spectra are not presented in the 
Appendix. 

Speed and elevation spectra from the same trial 
and site are placed together in the Appendix for comparison 
(Figures A-14 to A-23). It is evident from visual comparison 
that the -5/3 law fitted the rural speed and elevation 
spectra over a different range. Visual estimates were made 
of the low-frequency limit of agreement fF. between the 
spectra and a -5/3 power law. These estimates are given in 
Table 8. 

Examination of the table shows that the long- 
wavelength limit for rural speed spectra ranged between 10 
and 50 times the height. The long-wavelength limit for 
rural elevation spectra took on values between about 3 and 
6 times the height. The three rural spectra that agreed 
with the -5/3 law had long-wavelength limits of 8 to 12 


times the height. 
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It may be seen that the -5/3 law applied to one 
component when it did not apply to others. The range over 
which this occurred was from 1 to 4 octaves in the case of 
rural speed and elevation spectra. Local isotropy cannot 
occur when all three spectra do not obey the same power law. 
That is, agreement with a -5/3 law was found in the absence 
of local isotropy. Consequently, isotropy cannot be a 
necessary condition for a -5/3 law, and present theory must 


be considered incomplete. 
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CHAPTER 
CALCULATION AND COMPARISON OF TURBULENT INTENSITY 


5Ssde Methods) of Calculation 

Turbulent intensity is defined in terms of o/u, 
o/u, and oO /U, where 6, 0,, and o, are the standard 
deviations of the Cartesian velocity components u,v,w. 
Measures of turbulent intensity are also given by the 
standard deviation of the elevation angle, O43 the standard 
deviation of the azimuth angle, Oi and the ratio of the 
standard deviation of the speed to the mean speed, o./s- 
These quantities are measured in a spherical polar coordinate 
system. In this study the Cartesian turbulent intensities 
will be @pproxumatedsbyrithemmpolanm analogues-.gy Thatejis, 


it will be assumed that 


O/Saisamor/ us (52) 


These approximations are valid if the corresponding 
Cartesian and polar velocity components are related by the 


approximations used in section 4.2. 
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5.1.1 Choice of the Low-Frequency Limit for 
Variance Calculations 


Variances of azimuth, elevation, and speed were 
calculated from the low-frequency data. Exceptions were the 
azimuth variances for the trials on September 4 to 13, which 
were calculated from the high-frequency digitized data. 
in-thne calculation of variance from the low-frequency data, 
each record was divided into non-overlapping sections 656 
sec long. The variance was calculated from’each section 
using the mean calculated over the 656-sec period. The 
variances obtained from all the 656-sec sections were 
averaged. The sections were originally chosen to correspond 
to those used for power spectrum analysis. The length of 
the power spectrum analysis sections was later cut in half 
to improve the stability of the spectra. 

The calculated variances included energy in the 
spectrum above a frequency of approximately 1/656 Hz, or 
WW LX ge Hz. This frequency was. chosen to fall in the 
gap in the spectrum of horizontal wind speed observed by 
Van der Hoven (1957). His power spectrum covered the range 
£rom 2.x net Hz. to 0.25 Hz.) A gap in the spectrum was 
found between 5.6 x 105° HZpanG, 2/10: X iii HZesk, Gavan athe 
horizontal speed spectrum between 2.8 x LOMe HZ .ands5.16 x 
itp Hz was also observed by Walker (1964) under stable 
and neutral conditions. This gap may be interpreted as the 


gap between synoptic and micro-scales which must exist for 


a meaningful definition of atmospheric: turbulence. 
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The graphs of rural and urban power. spectra in the 
Appendix generally show decreasing power toward lower 
frequencies. This decrease would occur if a spectral gap 
was being approached. Exceptions to the decreasing trend 
occur in the rural azimuth and speed spectra which show 
Signs of low-frequency meandering, as discussed in sections 
4.2.2 and 4.2.3. These spectra pass through a minimum at a 
wavelength of about 200 m, and increase again toward lower 
frequencies. The lowest frequency resolved in the spectral 
analyses was 3.06 x soe HZ. 

On the basis of this observed trend of the power 
spectra, it may be concluded that most variances calculated 
in this study include essentially all the energy in the 
micrometeorological range. The exceptions are the azimuth 
and elevation variances for trials during which low-frequency 
meandering occurred. The variance may be underestimated 
in these cases. The urban azimuth variances for September 
4 to 13 must also be excepted. Only short sections of the 
time series were available for analysis in these cases. 
Consequently, the estimation of power at low frequencies 
was not stable. The contribution of low frequencies to these 
variances was estimated from the only reliable low-frequency 
urban azimuth spectrum available (September 20). 

5 leZ2en.Correction,ofalurbudent Intensity:,tor 

High-frequency Attenuation 

By examination of the spectra, it was found that, 


at wavelengths shorter than 4 m, the power decreased more 
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rapidly than predicted by the -5/3 law. This decrease was 
due to the limited frequency response of the Vector-Vanes 
discussed in section 4.1.2. At higher frequencies the 
spectrum was also attenuated by the low-pass filters. The 
calculated variances were corrected for this under estimation 
of the amount of energy at high frequencies by extrapolating 
the -5/3 law to all wavelengths shorter than 4 m. The 


energy above a frequency n, if the -5/3 law holds is given 


by 


D 


C -5/3 
Be 


Gry Oia ar 72 ne Sine) C52) 


°o 


The value of S(n) at 4 m which was used to make 
this correction was obtained by fitting a -5/3' line to the 
logarithmic spectra by eye and reading off the value of the 
ordinate at that wavelength. Because this value was in 
Gigltmzationiunits, it was multiplied by the calibration 
factor obtained by the method of section 2.4.1. The energy 
at wavelengths shorter than 4 m in the unmodified spectrum 
was removed by subtracting the sum of the power spectrum 
estimates above this wavelength from the variance, and the 
term 3/2 n,S(n,) was added. Table 9 contains the percentage 
corrections to the turbulent intensities that were made. 

To first order, the correction to the turbulent intensity 


is one-half of the percentage correction to the variance. 
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It may be seen from the table that the corrections 
applied to the turbulent intensities were all less than 10 
percent. “Corrections to the’ rural: vertical’ turbulent 
intensities were the largest. This is attributable to the 
larger proportion of energy at short wavelengths in the 
vertical spectrum. The smaller corrections to urban vertical 
turbulent intensities are explained by the shift of energy 
toward longer wavelengths in the city described in section 
A22.148 
Seles Correction of Urban Lateral Turbudent 

Intensities 

During four trials, the urban azimuth time series 
contained frequent discontinuities when the Vector-Vane 
azimuth potentiometer jumped from one end of its range to 
the other. Consequently, very little of the: record was 
available for analysis. Removal of the discontinuities 
from the digitized data was attempted. It was not successful 
because noise spikes, which were not removable by automatic 
means, frequently accompanied the discontinuities. Because 
of the volume of data to be processed, hand processing was 
not practical. The high-frequency digitized data were 
used for the calculation of variances in these cases because 
selection of the useable data was greatly simplified. The 
power spectrum and turbulent intensity were calculated from 
the low-frequency data for the September 20 trial because 


few discontinuities were present. 
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To achieve a reasonable level of significance in 
the calculated variances the long-wavelength energy was 
filtered out by calculating the variances over sections 41 
sec long. In this way energy below 1/41 Hz (0.024 Hz) was 
excluded. 

Numerical integration of the spectrum for the 
September 20 trial indicated that 30 percent of the energy 
was contained in wavelengths between 0.024 Hz and 0.00153 
Hz. A noise peak was evident in the high-frequency power 
spectra computed for September 4 to 13. The noise peak was 
responsible for about 7 percent of the variance. The 
high-frequency turbulent intensities for September 4 to 13 
were multiplied by 0.93 to compensate for the effect of the 
noise peak, and by 10/7 to compensate for the exclusion of 
low-frequency energy. The corresponding correction to the 
turbulent intensity was an increase of 15 percent. The major 
uncertainty in this correction is the amount of energy below 
0.024 4Hz. |}2E this amount. ranged from 20 percent to 40 
percent, the correction to the turbulent intensity would 
range “Lrom)S8 percent to 25 percent. Consequently, errors of 


more than 10 percent in the corrected values are unlikely. 


SieZ Comparison of Results 


The corrected values of turbulent intensity which 
were obtained by the technigues previously discussed are 
shown in Table 10. The urban turbulent intensity was 


greater than the rural turbulent intensity whenever concurrent 
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observations were available. The ratio of the average urban 
turbulent intensity to the average rural turbulent intensity 
was 4.6 for the vertical component, 3.7 for the longitudinal 


component, and 3.2 for thembaterakscomponent . 
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CHAPTER. VI 


CONCLUSIONS 


Turbulent intensities and power spectra at an 
urban and a rural site were measured in conditions that 
were generally stable with light winds. The scales and 
intensities of turbulence were compared and some systematic 
differences were found. Some aspects of the power spectra 
were compared, to the predictions. of theory. The conclusions 
which were drawn are summarized below. 

The intensity of turbulence was found to be 
more than three times as great at the urban site as at the 
mural. Site. 

A -5/3 power law fitted the urban power spectra of 
wind speed, elevation angle, and azimuth angle well at 
short wavelengths. The power law fitted these spectra at 
wavelengths shorter than 8 to 50 times the distance to the 
ground. 

Rural speed spectra were also well represented by 
a -5/3 law. The long-wavelength limit was 10 to 30 times 
the height above ground. Rural elevation spectra did not 
agree with a -5/3 law over as wide a range of wavelengths 
as did speed spectra or urban elevation spectra. The rural 
elevation spectra also exhibited more scatter in the region 


where agreement was obtained. The long-wavelength limit for 
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55 
rural elevation spectra was 3 to 8 times the height. Rural 
azimuth spectra agreed poorly with a -5/3 law. However, the 
rural azimuth signals were the smallest recorded, and it is 
likely that the signal-to-noise ratio was higher. 

A -5/3 law was observed to agree well with the 
shape of many spectra at wavelengths long in comparison 
to the height above the ground. Also, a -5/3 law applied to 
rural speed spectra over a range of 1 to 4 octaves when it 
did not apply to rural elevation spectra. It was concluded 
that isotropy is not a necessary condition in order that a 
-5/3 law should apply to atmospheric turbulence spectra. 

The peak of the elevation angle spectrum occurred 
at longer wavelengths at the urban site than at the rural 
Site. The peak of the rural elevation spectrum agreed with 
the near-neutral value of the universal function proposed 
for vertical velocity spectra by Busch and Panofsky (1968). 
This function was obtained on the basis of similarity theory 
and several sets of recent measurements at different 
locations. The peak of the urban elevation angle spectrum 
was an octave and a half lower than predicted by Busch and 
Panofsky. Two alternative hypotheses were proposed to 
explain the discrepancy. 

The first alternative is that vertical velocity 
spectra do not obey similarity theory close to the ground 
in a rough urban environment. That is, the vertical scale 
of turbulence in the city decreases more slowly with height 


than the linear relationship required by similarity theory. 
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It should be noted that observations in the city were made 
at a height small in comparison to the height of the major 
roughness elements. Consequently, an important assumption 
underlying the derivation of the similarity theory does not 
apply. 

The second. alternative is that the vertical velocity 
spectrum and the elevation angle spectrum are not similar in 
shape when the turbulent intensity is high. However, the 
agreement of the polar turbulence spectra with a -5/3 law 
suggests that the same relationships hold for polar and 
Cartesian spectra, at least at short wavelengths. 
Investigation of this problem is outside the scope of the 
present study. The literature provides little information 
on the relationship between turbulence statistics in polar 
coordinates and their Cartesian counterparts. Research on 
this topic would be of considerable interest. Until this is 
done, the first alternative hypothesis must be only 


tentative. 
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Figure A-l Elevation angle spectra, September 4 
(n*S(n) vs log £) 
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Figure A-2 
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Elevation angle spectra, September 5 
(neS(n) vs log f£) 
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Figure A-3 Elevation angle spectra, September 13 
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Figure A-4 Elevation angle spectra, September 20 
(n*8(n) ve log f) 
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Figure A-6 Rural speed and elevation spectra, September 17 
(neS(n) vs log f) 
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Figure A-7 Speed spectra, September 4 
(neS(n) vs log £) 


NxS (N} 


Figure A-8 


NxS (N) (X107!} 


68 


0.16 0.20 


le 


Q. 


0.08 


0.04 


.00 


"3.00 -2.50 -2.00 -1,50 -1,00 -0.50 0.00 0.50 
GG ar 
Urban 

oO 

w 

oS 


0.16 


0.12 


0.08 


ou 


,0.00 


3.00 -2.50 -2..00 -1.50 ~1.00 -0.50 0.00 0.50 
LOG F 


Rural 


Speed spectra, September 5 
(n°-S(n) vs log £) 
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Speed spectra, September 13 
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Figure A-9 
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Figure A-10 Speed spectra, September 20 
(nS (n'y va’ log £) 
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Figure A-1ll Rural azimuth spectra, September 4 and 5 
(n°S(n) vs idog £) 
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Figure A-12 Rural azimuth spectra, September 13 and 17 
(neS(n) vs log f) 
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Figure A-13 Urban and rural azimuth spectra, September 20 
(neS(n) vs log f) 
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Figure A-14 Urban speed and elevation spectra, September 4 
(log; és ve log" £) 
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Figure A-15 Urban speed and elevation spectra, September 5 
(log S vs log £f) 
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Figure A-16 Urban speed and elevation spectra, September 9 
(log S vs log £f) 
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Figure A-17 Urban speed and elevation spectra, September 13 
(log S vs log f) 
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Figure A-18 Urban speed and elevation spectra, September 20 
(log S vs log £f) 
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Figure A-19 Rural speed and elevation spectra, September 4 
(LOG S.°V.e. LOO. a) 
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Figure A-20 Rural speed and elevation spectra, September 5 
(log S vs log f) 
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Figure A-21 Rural speed and elevation spectra, September 13 
(log S vs log f) 
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Rural speed and elevation spectra, September 17 
(log S vs log £f) 
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Figure A-23 Rural speed and elevation spectra, September 20 
(log S vs log £f) 
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Figure A-24 Rural azimuth spectra, September 4 and 5 
(log S vs log f) 
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Figure A-25 Rural azimuth spectra, September 13 and 17 
(log S vs log £) 
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Figure A-26 Urban and rural azimuth spectra, September 20 
(log S vs log f) 


ii 1 
i”. 


a) 


Yaseen = 


; 1 2 
; te tree La) 


By is a ss Ba i ey ii ie " aA @ 
ty: ¥ a ioe J es a i ll Basi 
‘ ; f ‘ ely» - Ves ; 


tis eee 


4 a 
fe Oy, = 


- 7 . : 
"ve 
bert hil 
aie © 
= 


